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We discuss Faddeev-Popov quantization at the non-perturbative level and show that 
Gribov's prescription of cutting off the functional integral at the Gribov horizon does not 
change the Schwinger-Dyson equations, but rather resolves an ambiguity in the solution 
of these equations. We note that Gribov's prescription is not exact, and we therefore turn 
^ ■ to the method of stochastic quantization in its time-independent formulation, and recall 
d the proof that it is correct at the non-perturbative level. The non-perturbative Landau 

gauge is derived as a limiting case, and it is found that it yields the Faddeev-Popov 
method in Landau gauge with a cut-off at the Gribov horizon, plus a novel term that 
corrects for over-counting of Gribov copies inside the Gribov horizon. Non-perturbative but 
truncated coupled Schwinger-Dyson equations for the gluon and ghost propagators D{k) 
and G{k) in Landau gauge are solved asymptotically in the infrared region. The infrared 
critical exponents or anomalous dimensions, defined by D{k) ~ l/(/c^)^+"° and G{k) ~ 
1/(A;^)^^"^ are obtained in space-time dimensions d = 2,3,4. Two possible solutions 
are obtained with the values, in d = 4 dimensions, gg = 1, = —2, or ac = (93 — 
v^201)/98 ^ 0.595353, ao = -^ac- 



1. Introduction 

The problem of confinement in QCD presents a challenge to the theorist. One would 
like to understand how and why QCD describes a world of color-neutral hadrons with a 
mass gap, even though it appears perturbatively to be a theory of unconfincd and massless 
gluons and quarks. A basic insight into the origin of the mass gap in gluodynamics was 
provided by Feynman [1], Gribov [2], and Cutkosky [3]. These authors proposed that the 
mass gap is produced by the drastic reduction of the physical configuration space in non- 
Abelian gauge theory that results from the physical identification, Ai A2, of distinct 
but gauge-equivalent field configurations, A2 = ^Ai. A simple analogy is the change in the 
spectrum of a free particle moving on a line, when points on the line are identified mod- 
ulo 27r. The real line is reduced to the circle so the continuous spectum becomes discrete. 
In analytic calculations, the appropriate identification of gauge-equivalent configurations 
requires non-perturbative gauge-fixing, and in this article we approach the confinement 
problem by considering the how non-perturbative gauge-fixing impacts the Schwinger- 
Dyson equations of gluodynamics and their solution. We discuss both the Faddeev-Popov 
formulation, for which the Gribov problem may have an approximate - but not an exact 
- solution, and stochastic gauge fixing which overcomes this difficulty [4] . We also briefly 
compare our results for the infrared critical exponents with numerical evaluations of QCD 
propagators. In this connection we note that stochastic gauge fixing of the type considered 
here has been implemented numerically with good statistics on impressively large lattices 
(48^ X 64) [5] , [6] , [7] , and [8] . This opens the exciting perspective of close comparison of 
analytic and numerical calculations in this gauge. Stochastic quantization has also been 
adapted to Abelian projection [9], and convergence of the stochastic process has been 
studied theoretically [10]. 

We briefly review Faddeev-Popov quantization as a non-perturbative formulation. Wc 
note that the Faddeev-Popov weight Pfp(^) possesses nodal surfaces in A-space where 
the Faddccv-Popov determinent vanishes, and that a cut-off of the functional integral on a 
nodal surface does not alter the Schwinger-Dyson equations, because it does not introduce 
boundary terms. As a result, Gribov's prescription to cut off the functional integral at the 
(flrst) Gribov horizon [2], a nodal a surface that completely surrounds the origin [11], does 
not change the Schwinger-Dyson equations of Faddeev-Popov theory at all, but rather it 
resolves an ambiguity in the solution of these equations. 

We recall that Gribov's prescription is not in fact exact because there are Gribov 
copies inside the Gribov horizon that get over-counted. We then turn to the method of 
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stochastic quantization as described by a time-independent diffusion equation in ^-space 
(so tliere is no fictitious "fiftli time") [4]. Tliis metliod by-passes tlie Gribov problem 
of choosing a representative on each gauge orbit because gauge-fixing is replaced by the 
introduction of a "drift force" that is the harmless generator of a gauge transformation. 
We next derive a formulation of the Landau gauge, that is valid non-perturbatively, as a 
limiting case of stochastic quantization. It yields the Faddeev-Popov theory with a cut-off 
at the Gribov horizon, plus a novel term the corrects for over-counting inside the Gribov 
horizon. 

However attractive a formulation may be that is valid at the non-perturbative level, 
it would remain largely ornamental without actual non-perturbative calculations. Fortu- 
nately, progress in finding approximate but non-perturbative solutions for the propagators 
in QCD has been achieved recently within the framework of Faddeev-Popov theory both 
in Coulomb gauge using the hamiltonian formalism [12], and in Landau gauge by solving 
a truncated set of Schwinger-Dyson equations [13], [14], and [15]. The Schwinger-Dyson 
approach is reviewed in [16]. In the latter part of the present article we solve the Schwinger- 
Dyson equations in the non-perturbative Landau gauge to obtain the infrared critical ex- 
ponents or anomalous dimensions of the gluon and ghost propagators D{k) and G{k) in 
d = 2,3 and 4 space-time dimensions. The novel term is ignored here, in order to compare 
with other recent calculations, but we explicitly select the solution to the Schwinger-Dyson 
equations that vanishes outside the Gribov horizon. Although a truncation is necessarily 
required to solve these equations, nevertheless the values obtained for the infrared asymp- 
totic dimensions agree with exact results for probability distributions that vanish outside 
the Gribov horizon [17], namely the vanishing of D{k) at A; = 0, and an enhanced infrared 
singularity of G{k). These properties also characterize the non-perturbative solutions of 
the Schwinger-Dyson equations in QCD obtained in recent studies, [13], [14], and [15], 
and we verify that they have also adopted the solution that vanishes outside the Gribov 
horizon. 



2. Faddeev-Popov quantization at the non-perturbative level 

The standard Faddeev-Popov Euclidean weight in Landau gauge is given by 

(2.1) 
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with partition function 

Z{J)= f dAPFp{A) exp{J,A) 

\ (2.2) 

where (J, ^) = / cZ^xJ^A^. It depends only on the transverse components of J and we 
write Z = Z{J^'^). The complete set of Schwinger-Dyson (SD) equations reads 

^{S/SJ'^) Z = J^^ Z, (2.3) 

where "EiA^') = Sym{A^') - Tvln[-d^D^{A^')]. 

The Faddeev-Popov operator —d^D'^''{A^^') = —d^S°''^ — f°-^'^A^^'^d^ is hermitian be- 
cause A^^ is transverse, d^A^^ = 0. However it is not positive for every A^"^ . Because the 
Faddeev-Popov determinent is the product of non-trivial eigenvalues, det[—d^D ^^{A^^)] = 
Yin ^n{A^'^), it vanishes together with Q^piA^'^) whenever any eigenvalue vanishes, and 
the equation Xn{A^'^) = defines a nodal surface of Qfp(^*^) in A*''-space. The nodal 
surface where the lowest non-trivial eigenvalue vanishes, defined by Xo{A^'^) = 0, defines 
what is known as the (first) Gribov horizon. It forms the boundary of the region O, known 
as the Gribov region, with the defining property that all non-trivial eigenvalues of the 
Faddeev-Popov operator — 9^D^(A*'') are positive. Clearly Qfp(^*'^) is positive inside O 
and vanishes on dfl. It is known that the Gribov region has the following 3 properties: (i) 
it is a convex region of yl-space, (ii) it is bounded in every direction, and (iii) it includes 
the origin [11]. 

The existence of nodal surfaces of Qfp{A^'^) implies that the solution of the Schwinger- 
Dyson (SD) equation for Z{ J^'^) is not unique. For if the integral (2.2) that defines Z(J*'') 
is cut-off on any nodal surface, the same SD equation follows, without any boundary 
contribution. Moreover since the SD equation for Z{J^^) is linear, any linear combination 
of two solutions is a solution. These ambiguities are reflected in corresponding ambiguities 
in the solution of the SD equation for W{J^'^) = lnZ(J*'^), and for the effective action 
r(^*'^), obtained from W{J^'^) by Legendre transformation. 

[We illustrate these points by a baby model. Replace the field ^^^(^c) by a real variable 
a, and the Faddeev-Popov weight Qfp{A^^) by p(a) = exp(— |a^)(l— ^^a^), which vanishes 
on the "Gribov horizon" a = ±^~"^. The partition function Zi(j) = da p{a) ex.p{ja) 
satisfies the SD equation 

id) [|-(^-4)"'^^(-4)"1 '^-^^ 



corresponding to the action S = — ln(l — g^o?). Suppose we restrict the integral 
to the Gribov region, \a\ < g~^, so the partition function is given instead by Z2{j) = 
f^y^da p{a) exp{ja). The change occurs only for > so Zi{j) and Z2{j) have 

the same perturbative expansion. It is clear that Zi{j) and Z2{j) satisfy the same SD 
equation (2.4) without a boundary contribution, because p{a) vanishes on the boundary 
a = ±1/^. Moreover, because the SD equation for Z{j) is linear, any linear combination, 
Z{j) = aZi{j) + pZ2{j) also satisfies the same SD equation. Of course only one of them 
corresponds to a weight that is positive everywhere. This example easily generalizes to 
any number of dimensions.] 

Gribov proposed to cut off the integral on the boundary dfl of what we now call the 
Gribov region Q, so the partition function is given by 

Zn{J'')= I c^^^■^QFp(^'Oexp(J*^^'■^), (2.5) 
Jq. 

for he conjectured that the region where the Faddeev-Popov operator is positive contains 
only one gauge copy on each gauge orbit. Since dVt is a nodal surface, satisfies the same 
Schwinger-Dyson equation (2.3), j^{5/5J^'^)Zq, = J^'^Zq,. Instead, Gribov's proposal 
selects a particular one out of a class of non-perturbative solutions of these equations.^ 

It is now known however that Gribov's conjecture is not exact. Indeed, there are 
Gribov copies inside the Gribov horizon.^ We shall show however that an exact non- 
perturbative formulation yields Gribov's proposal plus a well-defined correction term that 
corrects for overcounting inside the Gribov horizon. 

^ For partition function Zn, the expectation value of J^^ in the presence of the source J**^ 
is Ac\ = {A^^) ju = 6Wn/SJ^'^ , where Wn = InZn. Because the probability distribution in 
the presence of sources (5(A*'^) exp( J*"^, A**^) is positive in the Gribov region and because 
is convex, it follows that Ad lies in fi. Consequently the effective action rn(Aci), obtained by 
Legendre transform from Wn, is defined only on the Gribov region. 

^ As shown in sec. 4, they are given by ^A, where g = fifmin(a^) is any local minimum of the 
functional FA{g) = J d'^x \^A\'^, where ^A = g^^Ag + g^^dg is the transform of A by the local 
gauge transformation g{x). In a lattice discretization, the link variable corresponding to the 
field A{x) is generically a random field, so the minimization problem is of spin-glass type which 
is known to have many solutions. On the other hand for a smooth configuration, such as the 
vacuum, A = 0, there are few solutions. Thus the number of copies is different for different orbits. 
Moreover, since the Faddeev-Popov weight is positive inside the Gribov horizon, there can be 
no cancellations to save the day. Note also that, in a lattice discretization, the variables that 
characterize a configuration take values in a compact space, so a minimizing configuration always 
exists, which shows that Q contains at least one Gribov copy for each orbit. 
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There is an alternative proposal to make the Faddeev-Popov method valid non- 
perturbatively. It is conjectured that if one sums over all Gribov copies using the signed 
Faddeev-Popov determinent det[— then additional Gribov copies cancel in 
pairs, the reason being that the signed determinent counts the signed intersection number 
which is a topological invariant. This is presumably the outcome of BRST quantization 
which has formal properties that suggest it may be valid non-perturbatively [18]. Moreover 
this conjecture is supported by simple models [19] and [20]. However it is not known at 
present how to turn this prescription into a non-perturbative calculational scheme in QCD, 
for example, by selecting a particular solution of the SD equations, and moreover if the 
measure is not everywhere positive, there is the danger of delicate cancellations that may 
cause an approximate solution to be unreliable. On the other hand, the Gribov proposal is 
easily implemented, for example by requiring that the solution of the SD equation possess 
positivity properties. 

3. Time-independent stochastic quantization 

The difficulties with Faddeev-Popov gauge fixing pointed out by Gribov are by-passed 
by stochastic or bulk quantization. This method is a formalization of a Monte Carlo sim- 
ulation [21], and in its most powerful formulation it makes use of a fictitious "fifth time" 
that corresponds to computer time or number of sweeps of the lattice in a Monte Carlo 
simulation. Despite the Gribov ambiguity, there is no problem of over-counting with gauge 
fixing in Monte Carlo simulations, which is achieved by a gauge transformation of choice 
after any sweep, nor is there one in stochastic or bulk quantization which relies on an in- 
finitesimal gauge transformation. In the 5-dimensional formulation, the field = t) 
depends on the 4 Euclidean coordinates and on the fifth time t. One may easily write 
down the SD equations for this formulation which involves a local 5-dimensional action, 
and BRST invariances and Slavnov- Taylor identities are available to control divergences 
[22], [23], and [24]. However the 5-dimensional propagator D = D{k'^,u) depends on two 
invariants k'^ and cu, which makes the solution of the SD equations in 5 dimensions more 
complicated, and we shall not attempt it here. 

Instead we turn to an older 4-dimensional formulation of stochastic quantization [4]. 
It is based on an analogy between the (formal) Euclidean weight Pym(^) = exp[— ^ymI^)] 
and the Boltzmann distribution P{x) = cxp[— F(a;)]. The latter is the solution of the time- 
independent diflFusion equation + §^)P = 0) where the drift force is Ki = 
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We shall shortly consider more general drift forces Ki that are not necessarily conservative. 
The field theoretic analog of this equation is 

MX).-/.S^(^ + ^)P(X) = 0. (3.1) 

where the drift force is KYM,ij.ix) = — ff^^^ ; which is solved by P{A) = exp(— S'ym)- 

For a gauge theory, this solution is not normalizable. However for a gauge theory, 
one may modify the drift force Kym Kym + Kgt by adding to it a "force" Kgt tangent 
to the gauge orbit, without changing the expectation- value of gauge-invariant observables 
0{A). Such a force has the form of an infinitesimal gauge transformation -f^gt,^ = {D^vY^ 
where v°'{x\A) is an element of the Lie algebra, and {D,j_vY = d^v"" + /'^''^A^^v^ is its 
gauge-covariant derivative. This force is not conservative, which means that it cannot be 
expressed as a gradient, -K'gt.p 7^ ~jx~-> ^^^^ method is not available in a local action 
formalism in 4-dimensions. The total drift force is given by 

-f^M = KYM,fj, + 

and -P(A) is the solution of the modified time-independent diffusion equation 
EP={HYM + ii^t)P 

It is easy to show [4] that the expectation- value (O) = f dA 0{A)P{A), of gauge-invariant 
observables 0{A) is independent of v, using the fact that H^^ is the generator of an 
infinitesimal local gauge transformation 

= G{v) = I dt^x v{x)G{x) 

and that 0{A) and Hym are gauge invariant, G{x)0 — 0, and [G{x), Hym] = 0. 

The additional drift force Kgt^/j, = -D^^v must be chosen so that it is globally a restoring 
force along gauge orbits, thus preventing the escape of probability to infinity along the 
gauge orbit where Sym is flat. This may be achieved by choosing Kg^ to be in the 
direction of steepest descent, restricted to gauge orbit dirctions, of some conveniently 
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chosen minimizing functional F(A). A convenient choice is the Hilbert square norm, 
F{A) = ll^lp = J d'^x\A^\'^. For a generic infinitesimal variation restricted to gauge 
orbit directions SA^ = eD^uj, we have 

5F{A) = 2{A^,5A^) = 2{A^,eD^uj) = 2e{A^,d^u;) = -2e{d^A^,oj). (3.5) 

The direction of steepest descent of restricted to gauge orbit directions, is seen to 

be SAf^ = eD/j^u, for u; = dxAx. Thus if we choose v = a~^dxAx, where a is a positive 
gauge parameter, the drift force -f^gt,^^ = ^"^-D^^a^a points globally in the direction of 
steepest descent, restricted to gauge orbit directions, of the minimizing functional 
In the following we shall use the time-independent difi^usion equation 

The 5-dimensional formulation is based on the corresponding time-dependent diffusion 
equation 

dP/dt = -HP. (3.7) 

4. Non-perturbative Landau gauge 

Because the gauge-fixing force points in the direction of steepest descent of the mini- 
mizing functional F{A) = 1 1^| |^, restricted to gauge orbit directions, it follows that for large 
values of the gauge parameter the probability gets concentrated near the local minima 
of this functional restricted to gauge orbit variations.^ At a local minimum the first varia- 
tion vanishes for all 5F{A) = — 2e(a;, d^A^) = 0, as we have just seen, so at a minimum, 
the Landau gauge condition d^A^ = 0, is satisfied. In addition, the second variation in 
gauge orbit directions is non- negative, 5'^F{A) = — 2e(a;,9^M^) = — 2e^(a;, 5^£)^a;) > 0, 
for all u, which is the statement that the Faddeev- Popov operator M{A) = —dfj,Df^{A) is 
positive. These are the defining properties of the Gribov region, and we conclude that in 
the limit in which the gauge parameter approaches zero, o ^ 0, the probability P{A) gets 
concentrated on transverse configurations A = A^'^ that lie inside the Gribov horizon. We 
have noted above that there are Gribov copies inside the Gribov horizon. However the 
present method does not require that the probability gets concentrated on any particular 

^ These conditions define a local minimum at g{x) = 1 of the functional on the gauge orbit 
through A defined by Fa (5) = j d'^x\<^A\^. 
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one of them such as, for example, the absolute minimum of the minimizing function, and 
for finite gauge parameter a, the gauge- fixing is "soft" in the sense that no particular gauge 
condition is imposed. For gauge-invariant observables, it does not matter how the proba- 
bility is distributed along a gauge orbit, but only that it be correctly distributed between 
gauge orbits. This is assured because a harmless gauge transformation was introduced 
instead of gauge fixing. 

We have noted that A becomes purely transverse in the limit a — > 0. We shall solve 
(3.6) in this limit by the Born-Oppenheimer method in order to obtain the non-perturbative 
Landau gauge. For small a, the longitudinal component of A is small and, as we shall see, 
it evolves rapidly compared to the transverse component. However because of the factor 

in (3.6), the mean value of the longitudinal part of the gluon propagator strongly 
infiuences the transverse propagator in the limit a ^ 0. 

We decompose A into its transverse and longitudinal parts according to A^^ — A^^'^ + 
al/2a^(a2)-lL^ so a^Aj = al/2L^ and ^ = - a'^l'^d^-^. In terms of these 

variables, eq. (3.6) reads 



where we have used the notation {K x L)^ = [K, L]^ — f'^^'^K^L'^, for elements K and L 
of the Lie algebra. The leading terms in H are of order a~^, a~^/^ and aP , and we expand 
P = Pq + a^/^Pi -|- ... . The leading term, of order a~^, reads 

j d'x ji[{-d^)ji - d,D,{A'lL) Po = 0. (4.2) 

This is solved by Pq that is Gaussian in L, 

PoiA'\L) = Q{A'') (detX)-^/2 exp[-l/2(L, XL)], (4.3) 

where X^'^(a;, y; A*^) is a symmetric kernel. Equation (4.2) is satisfied provided X satisfies 
{L,X{-d'^)XL) - {L,XML) = identically for all L, and tr[{-d'^)X - M] = 0. Here 
M = M{A^^) = —d^D^{A^^) is the Faddeev- Popov operator that is symmetric for A 
transverse. The first equation yields X{-2d'^)X = XM + MX, or MY + YM = -28'^ for 
Y = X~^. Moreover when this equation is satisfied, it implies that the second equation 
is also satisfied. The equation for Y is linear. To solve it we take matrix elements in the 
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basis provided by the eigenfunctions of the Faddeev- Popov operator M{A^^)un = A„tt„, 
where An = Xn{A^'^), and obtain 

{Um, X~^Un) = {Un, YUm) = {Xm + K)~^{Un, {-2d'^)Um)- (4.4) 

We see that the Gaussian solution Po{A^^, L) is normahzable in L only when all the eigen- 
values XniA^"^) are positive, namely, for A^^ inside the Gribov region. However we have 
seen above that in the limit a — > the solution P{A) is supported inside the Gribov region. 
Thus the coefficient function Q{A^'^) carries a factor ^(Ao(^*^'^)), that restricts the support 
of Pq to this region. Finally we note that for A^"^ in the Gribov region, Y may be written 

j-OO 

X-^=Y= dt exp{-Mt){-2d'^)exp{-Mt). (4.5) 
Jo 

This representation shows explicitly that X is a positive operator for A^^ inside the Gribov 
region. 

To determine Q{A^^) we substitute (4.3) into (4.1), and integrate over L. This kills 
the term in j^. It also kills the term of order a~^/^ in a^^/^A^ x L = a^^/'^[A^^ + 
a^^'^d^{d'^)~^ L] X L because this term is odd in L. This gives in the limit a — > 0, the finite 
equation for Q{A^'^), 

Here -K^gteff is the average over L of the gauge-transformation force, with weight 
(detX)-V2 exp[-l/2(L,XL)], namely 



= f^%{dY'y'"'{x,y;A'^)\y=, (4.7) 



L 



dtf''%{dY^[eM-Mt){-2d^)eM-Mt)]''%,y)\ 



y=x-, 



and we have used Y'^'^{x^y; A^^) = {L^{x)L^{y)). We now take the limit a — > namely 
P{A) lima^o P{A) = lima^o ^'o(^)- With L = a-'^/'^d^.A^ this gives 

P{A) = Q{A'^)S{d^A^), (4.8) 



where Q{A^^) is the solution of (4.6). This defines the non-perturbative Landau gauge. 
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To exhibit the relation between the non-perturbative Landau gauge and Faddeev- 
Popov theory, we decompose i^^gteff, given in eq. (4.7), according to 

i^gteff ^Ki+K2 

l-OO 

Kl,{x) = / dt f^%id') 
Jo 

The first term may be written 

so Ki is a conservative drift force, derived from an action S = tr In M = — (Indet M) that 
precisely reproduces the Faddeev-Popov determinant. So if K2 were neglected we regain 
the Faddeev-Popov theory, with the added stipulation to choose the solution that vanishes 
outside the Gribov horizon. 

The second term may be simplified using the identity 

[a^exp(-Mt)] = - /" ds exp{-Ms)[d'^,M]exp[-M{t- s)], (4.11) 
Jo 

which gives 

l>00 

Klj^) = - / ds f'^%{dY'{eM-Ms)[d^M] eM-Ms)M-'r^{x,y)\y^,, 
Jo 

(4.12) 

where, M = M(A*''). The "drift force" K2 is a novel term. It's presence is required to 
correct the overcounting, discussed in sec. 2, that occurs when the Faddeev-Popov theory 
is cut off at the Gribov horizon."^ 



-^{[2d'',eM-Mt)]eM-Mt)r''{x,y)\y=,. 



5(trlnM) 



tr,b/ 



(4.10) 



5. Schwinger-Dyson equations 

The partition function is defined by 

Z{J) = J dAQ{A"^)S{d^A^)eMJ,A) 



r (5.1) 

= / dA^' Q{A'')e^^{J,A^'). 



^ In our derivation we used the Born-Oppenheimer method that is non-perturbative in g in 
order to obtain the a — >^ limit at finite g. So the presence of the new term K2 is not in con- 
tradiction with the fact that the Faddeev-Popov theory provides a formal perturbative expansion 
that has all the correct properties including perturbative unitarity. 
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It depends only on the transverse component of (on-shell gauge condition), and we 
write Z = Z{J^^). Generally, in the Faddeev-Popov approach, one relaxes the transversal- 
ity condition, by writing 5(5^^^) = J dbexp{i J df^x h 9^^^), and then one uses Slavnov- 
Taylor identities to determine longitudinal parts of vertices. However these identities have 
not yet been derived in the present 4-dimensional stochastic approach, and we shall solve 
the SD equations using the on-shell formalism for the gauge-condition. The on-shell cor- 
relation functions, such as propagators, are the same as the off-shell ones, but the vertices 
(one-particle irreducible functions) are strictly transverse. Renormalization theory is not 
well articulated at present in the on-shell formalism, but we shall not encounter ultraviolet 
divergences in the SD equations in the infrared limit. Moreover, we shall see that in this 
limit the SD equations are invariant under the renormalization-group.^ 

The partition function Z(J*'^), which is the generating functional of (transverse) cor- 
relation functions, is the fourier transform of the probability distribution Q{A^^). Con- 
sequently the SD equation for Z(J*^) is simply the diffusion equation (4.6), expressed in 
terms of the fourier-transformed variables, 

j d'x J'; [j^^ + Ktoteff.M (^)] Z{Jn = 0. (5.2) 
Here we have introduced the total effective drift force 

Ktote^AA'l ^ -^-^(A'l + K,tef^AA'l. (5.3) 

Only the transverse component of -fCtotefF,/:^ appears in the following. The free energy 
W{J^'^) = lnZ(J*^), which is the generating functional of connected correlation functions, 
satisfies the SD equation 



/ 



d'^x J*; 



The effective action is obtained by Legendre transformation, r{A^^) = {J^^,A^^) — W{J^^) 
by inverting A^'^ = It satisfies the SD equation 



/ 



4 5T{A'')lST{A'') 



SAi 



SAi 



+ i^toteff ,M {a'' + T^i^'') sJ^)]=^^ (5-5) 



^ It should also be noted that it is not known at present how to maintain the Slavnov- Taylor 
identities exactly at the non-perturbative level in the off-shell formalism, although methods for 
dealing with this have been proposed [25]. 
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where the argument of i^^totefF is written in matrix notation, and is given exphcitly by 

1:^. Here 'Dni,{x,y; A^'^) — jjp^ 



+ J d^y Vfj,ty{x,y;A^'')jjiT-^. Here Vfj,„{x,y; A^'') = I^^fM is the gluon propagator 



in the presence of sources. 

To obtain the SD equation for the propagator, we expand in powers of A^'^ , 

"^^/^ (5.6) 

5 



5A^' 

where we have again used matrix notation. Here D^j^ = D^i,{x — y) is the gluon propagator 
in the absence of sources, and 

Both D and R are identically transverse, and in momentum space, by virtue of Lorentz 
invariance, are of the form D^^{k) = D{k'^)[6^i, — k^k^j/k"^] and -R^j^(/c) = — 
kjj.kiy/k'^]. Upon equating terms quadratic in A^^ in (5.5) we obtain 

J d'x [D-'A'^] (x) [{D-' - R)A';] (x) = 0, (5.8) 

which we write in matrix notation as 

{D-^A'\ [D-^ - R]A'') = 0. (5.9) 

This holds identically in A^^. Prom the expressions for D^v{k) and R^j^{k), we see that both 
are symmetric operators that commute, DR = RD. As a result, the operator appearing 
in the last equation is symmetric and must vanish, 

D-\D-^ - R) = 0. (5.10) 

This gives the SD equation for the gluon propagator 

D-^ = R, (5.11) 

where R is given in (5.7). 
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6. Solution of SD equation in the infrared 

Recall the decomposition -ftTgteff = Ki + K2, where Ki = is the drift force 

that, in the absence of K2, describes the Faddeev-Popov theory in Landau gauge. Since it 
is not without interest to solve Faddeev-Popov theory non-perturbatively in Landau gauge, 
and in order to compare our results with other authors, we shall here ignore K2, the novel 
term that corrects the over-counting that occurs when the Faddeev-Popov theory is cut 
off at the Gribov horizon. 

The remaining drift force, Ki, describes Faddeev-Popov theory in Landau gauge. We 
have seen in sec. 2 that there is an ambiguity in the solution of the SD equations of 
the Faddeev-Popov theory, with no clear prescription to resolve it at the non-perturbative 
level. Fortunately the present derivation provides the additional information that is needed 
to resolve this ambiguity: we must choose the solution of the SD equations that vanishes 
outside the Gribov horizon because, as we have seen, Q{A^^) vanishes outside the Gribov 
horizon in the limit a — > 0. With this choice it is likely that qualitative features of the 
exact theory (with K2) will be preserved. 

With neglect of K2 we may write directly the familiar SD equations of the Faddeev- 
Popov theory in Landau gauge, in an arbitrary number d of Euclidean dimensions 



where G{p) is the ghost propagator, ri,{p, k) is the full ghost-ghost-gluon vertex. In Landau 
gauge, a factorization of the external ghost momentum occurs, so the ghost-ghost-gluon 
vertex is of the form r,^{p, k) = r^^\{p, k)p\. As a result there is no independent renor- 
malization of (p, k) , and the renormalization constants in Landau gauge are related by 
Z^Z^Z"^ = 1, where ^fo = ZgQr, Dq — Z^Dr, and Go = Z^Gr- So far, we have written the 
SD equations for unrenormalized quantities, with the index suppressed. 

We must select the solution to these equations that corresponds to a probability 
distribution Q{A^'^) that vanishes outside the Gribov horizon. To do so, it is sufficient to 
impose any property that holds for this distribution, provided only that it determines a 
unique solution of the SD equations. Besides positivity, which will be discussed in the 
concluding section, there are two exact properties that hold for a probability distribution 




O G{p + k)ip + k)^ G{p) r,,{p, k) + (gluon loops) 



(6.1) 




d'^k G{k+p) {k + p)^D^,{k) r,ip,k), 



(6.2) 
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P{A^'^) that vanishes outside the Gribov horizon: (i) the horizon condition, and (ii) the 
vanishing of the gluon propagator at A; = 0, [17]^ 



limD(/c) = 0. (6.3) 

The horizon condition (i) is equivalent to the statement that G(j)) diverges more rapidly 
than or 

limyCip)]-^ =0. (6.4) 
Indeed if we divide the SD equation (6.2) by p^, and impose this conditon, we obtain 

6^x = Ng\27r)-''J d'^k G{k) D^,{k) V,xM. (6.5) 

This is the non-perturbative statement that the ghost self energy, which is of the form 
S(p) = p^E^a(p)pa because of the factorization of the external ghost momentum, exactly 
cancels the tree level term at p = 0, 

<^M = 5^ma(0). (6.6) 

Equations (6.5) and (6.6) are the form of the horizon condition given in [26], [27], and [28].^ 
We will see that it is sufficient to apply either condition (i) or (ii), and the other condition 
then follows automatically. The horizon condition allows us to write the SD equation for 
the ghost propagator, (6.2), in the form 

G-\p) = Ng\27r)-'' J d'^k p^D^,{k) [T,,x{0,k) G{k) - r,,x{p, k) G{k + p)]px, (6.7) 

where we have used k^D^^{k) = 0. This equation was solved numerically in 3-dimenions 
in [28], using an assumed form for D{k). 

We wish to determine the asymptotic form of the propators at low momentum, 
G^{p'^), and D^l{k) = D'^{k'^)P^^^{k), where Pj^Hk) = 5fj,„ - k^k„/k'^ is the transverse 



^ The vanishing of the gluon propagator at A; = results from the proximity of the Gribov 
horizon in infrared directions. 

^ In a space of high dimension N the probability distribution within a smooth surface such as 
a sphere r < R gets concentrated near the surface r = R because of the entropy or phase-space 
factor r^~^dr. The horizon condition is the statement that the probability distribution within 
the Gribov horizon is concentrated on the Gribov horizon because the dimension N of A-space 
diverges with the volume V in, say, a lattice discretization. 
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projector. For this purpose we let the external momenta in the SD equations be asymp- 
totically small compared to QCD mass scales. In this case the loop integration will be 
dominated by asymptotically small loop momenta, so the propagators inside the integrals 
may also be replaced by their asymptotic values. This is true provided that the result- 
ing integrals converge, as will be verified. We shall also truncate the SD equations by 
neglecting transverse vertex corrections, as usual, in order to obtain a closed system of 
equations, T^J{p, k) —>■ P^^^{k)p^. Such truncations may, possibly, be justified a posteriori 
by calculating corrections to see if they are small. Because Dn^{k) is transverse, the SD 
equation for the ghost propagator simplifies to 

J (6.8) 

X [ G^{k'^) - G^{{k^pf)]. 

This equation is invariant under renormalization because of the identity Z^Z^Z^ = 1. 
This allows us to take all quantities in the last equation to be renormalized ones, with 
suppression of the index r. 

Because the asymptotic infrared limit is a critical limit, the asymptotic propagators 
obey simple power laws, 

according to standard renormalization-group arguments. Here and ao are infrared 
critical exponents or anomalous dimensions that we shall determine, while |U is a mass 
scale, and cd and cg are dimensionless parameters. The horizon condition (6.4) implies 
aa > 0, whereas (6.3) implies ao < —1- Upon changing integration variable according to 
kjj, = \p\k'^, and equating like powers of p we obtain 

aD + 2aG = - {4:-d)/2. (6.10) 

The integral is ultraviolet convergent for (i— 2(1+q;o) — 2(H-aG) — 2 < 0, where 2 powers of 
A; are gained because of the difference [G^^{k'^)—G^^{{k+p)'^)]. Withao = — 2q!g — (4— d)/2, 
this gives ckg < 1 as the condition for ultraviolet convergence, so < ckg < 1. 

We now turn to the SD equation for the gluon propagator (6.1). In the exact Faddeev- 
Popov theory with off-shell gauge condition, the right hand side of (6.1) is exactly trans- 
verse in k on both free Lorentz indices /i and v by virtue of the Slavnov- Taylor identities. 
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This allows us to apply transverse projectors P^^i^{k) = S^u—kf^ku/k"^ and Pl^,^{k) to these 
indices. In our derivation, with on-shell gauge condition, the projectors are automatically 
applied. As a result, since the gluon propagators are transverse, only the transverse parts 
of the vertices contribute on the right-hand side. We therefore make the truncation ap- 
proximation of replacing these transverse vertices by their tree-level expressions. We now 
estimate the various terms on the right hand side of the SD equation (6.1) for D{k). We 
just concluded from the horizon condition and the SD equation for G{p) that aa > and 
a/) < 0. As a result, on the right-hand side of (6.1), the ghost loop that we have written 
explicitly is more singular in the infrared than the gluon loops. Morover in the infrared, 
D~^{k) ~ (A;^)*^^+"-°) is more singular at /c = than the tree-level term ~ k"^ because 
CKD < 0. We now let the external momenutm k have an asymptotically small value, so the 
loop integration is dominated by asymptotically small values of the integration variable p 
(provided the resulting integral converges). We take the asymptotic infrared limit of (6.1) 
with external projectors and obtain 

[D-'{k^)]-'Pll{k) = Ng\27T)-%\{k) J d'^p p^ G-%{k +pf) G-^p") p, P'^Hk). 

(6.11) 

We take the trace on Lorentz indices and obtain 



[D^^{e)]-^ = Ng^{2Ti)-\d - l)/c']-' j d'^p \p^k^ - {p ■ kf] ^^^((/c + pf) G'^'ip^). 

(6.12) 

Like the ghost equation (6.8), this equation is invariant under renormalization because of 



the identity ZgZ^Z^ = 1, and we may again take all quantities to be renormalized, with 
suppression of the index r. We substitute the power-laws (6.10) into this equation. By 
the power counting argument that was used for the ghost propagator, we again obtain the 
relation of the infrared critical exponents an + 2aG = — (4 — d)/2. This integral converges 
in the unltraviolet for d — 2 < 4(1 -|- ckg), or ao > {d — 2)/4. 
The gluon and ghost SD equations now read 

(cd4)-' = /d(«g) = /g(«g), (6.13) 

where 

loiao) = Ni27T)-^id - l)-i(fc2)-(2+«-) 

f (6.14) 

X / d'^p \p^k^ - {p ■ kf] [{k + p)2]-(1+«g) [^^)2p(l+aG) 
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J (6.15) 

and it is understood that ckd = — 2q;g — (4 — ci)/2. The critical exponent ao is determined 
by the equahty (6.13). The integrals Id{olg) and /g(q;g) are evaluated in the Appendix, 
without angular approximation, in arbitrary Euclidean dimension d. 



7. Determination of infrared critical exponents 

To determine the critical exponent ao, we substitute the formulas for lD{aG) and 
Iq^cx-g)-, given in the Appendix, into the equation loisx-G) = Ig{(^g) and obtain, for 
a = ckg, 

(rf-l)7r r(l + 2a) T{d-2a) _^ 

" sin(7ra) r(-2a + d/2) r(l + a + d/2) V{-a + d/2) r(l + 2a - d/2) ~ 

(7.1) 

We take the dimension d of space-time in the interval 2 < d < A. The integrals Id{ol) and 
Id{(x) are both convergent in the ultraviolet only for a in the interval < ((i— 2)/4 < a < 1, 
so the equation which determines a holds only in this interval. However whereas Id{c() 
is manifestly positive throughout this interval, the expression for /g(ci) is negative for 
a > d/4, because l/r(— 2q: + d/2), changes sign at o: = d/A, so we look for a solution only 
in the reduced interval 

< (d-2)/4< q; < d/4< 1. (7.2) 

The identity, r{-2a + ci/2)r(l + 2a- d/2) = n/sm[n{-2a + d/2)], gives 

(d-l)sin[7r(-2« + d/2)] r(l + 2a) r{d-2a) ^ 
J^^"^^- sin(7rQ;) T{l + a + d/2) T{-a + d/2) " ^ ^ 

For the case of physical interest, d = 4, the allowed interval is 1/2 < a < 1, and the 
function /4(a) contains the factor sin[7r(— 2q; + 2)]/ sin(7ra;) = sin[7r(— 2q;)]/ sin(7rQ;) which 
is of the indeterminate form 0/0 at a = 1. To control this (and a similar indeterminacy for 
d = 2 at a = 0), we first consider d in the range 2 < d < 4, and then take the limit d ^ 4 
(and d ^ 2). For d in this range, one sees that the function /^(q;) is positive and finite, 
fd{ci) > 0, for a in the interior of the allowed interval {d — 2)/4 < a < d/4, but vanishes 
at both end-points fd[{d — 2)/4] = /d(d/4) = because of the factor sin[7r(— 2q; + d/2)]. It 

17 



follows that the equation fd{a) = 1 has an even number of solutions (if any)^ for 2 < d < 4. 
We now set d = 4 and obtain, 



-3sin(27ra) r(l + 2a) r(4 - 2a) _ ^ 
"'^^"^ " sin(7ra) r(3 + a) r(2 - a) ~ 



(7.4) 



(7.5) 



We use 

r(l + 2a) r(4 - 2a) = (3 - 2a)(2 - 2a){l - 2a) 2a n/ sin(27ra) 
r(3 + a) r(2 -a) = {2 + a){l + a)a (1 - a)7r/ sin(7ra), 
and obtain 

(3-2a)(2«-l) 
= (2 + a)(l + a) = ^' (^-^^ 

where we have used = 1, which is valid only for 1/2 < a < 1. This yields 

a quadratic equation with roots a = [93 ± y(1201)]/98 [93 ± 34.66]/98. Only one root 
a 0.5953 lies in the interval 1/2 < a < 1. On the other hand we have just seen that for 
2 < d < 4, there are an even number of roots. The resolution is that for d = 4 — e there 
are two roots, and the second root is given by a = 1 — 0(e), so in the limit d ^ 4, there 
is a second root at o; = 1. 

We conclude that the infrared critical exponents a = aQ and «£> = — 2q;£) — (4 — (i)/2 
are given, in d = 4 dimensions, by two possible sets of values 

ac = 1; OLD = ~2 

ao = [93 - y(1201)]/98 ^ 0.5953 ao = -[93 - ^(1201)1/49 ^ -1.1906. ^^'^"^ 

In the same way one finds for cZ = 2, 

ao = 0; aD = —1 

(7.8) 

ao = 1/5; an = —7/5. 

For (i = 3, one obtains the equation, /^{a) = ^^"^(■3^2qx"i+2q)^""'^ ~ ''" ^^^^ roots in the 
interval 1/4 < a < 3/4, given by 

an = 1/2; ckd = —3/2 

(7.9) 

ac- 0.3976; ao ^ -1-2952. 



^ Prom numerical plots it appears that for 2 < d < 4 there are always 2 distinct real roots in 
the range (d — 2)/4 < a < d/4, except possibly 2.662 where there may be a double root 

near a ^ 0.33095. 
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We expect that in each case one of the roots is spurious, and arises because (7.3) does not 
express the full content of the theory. 

We note that in each case, one solution corresponds to ao = {d — 2)/2 = 0, 1/2, and 
1, for d = 2,3, and 4, which gives G{k) ~ l/(/c^)^/^. This may be too infrared singular 
to be acceptable. But for d = 2, the other solution, with ac — 1/5, is even more infrared 
singular, which suggests that for c? = 2 the first solution may be preferred namely, = —1 
and ao = 0, which may make the case d = 2 pathological in the Landau gauge. This case 
is exactly solvable in the axial gauge because the non-linear term in the d = 2 Yang-Mills 
field is absent in this gauge, and gives an area law at the classical level. There can of 
course be no physical gluons in d = 1 + 1 dimensions even in the free theory which may 
thus be considered confining. Clearly the case d = 2 in the Landau gauge requires a more 
detailed investigation that we do not attempt here. 

8. Discussion and Conclusion 

We have seen that because the Faddeev- Popov weight Pfp{A) contains nodal surfaces, 
the SD equations corresponding to the Faddeev-Popov method are ambiguous, and in 
practice one does not know how to select an exact and globally correct solution. Gribov's 
proposal, to cut-off the Faddeev-Popov integral at the first nodal surface, produces a 
positive probability distribution, but it is not exact because it overcounts some gauge 
orbits, although it may give a useful approximation. 

By contrast the method of stochastic quantization by-passes the Gribov problem of 
selecting a single a representative in each gauge orbit. Instead the diffusion equation 
in A-space, eq. (3.6), contains an additional "drift force" a~^D^ d ■ A that is a harmless 
generator of a gauge transformation. The corresponding DS equation that defines the non- 
perturbative Landau gauge was obtained by solving the limit a ^ of this equation by the 
Born-Oppenheimer method. The limiting probability distribution Q{A^^) was shown to 
vanish outside the Gribov horizon. It is determined by a diffusion equation that contains 
the novel term K2, eq. (4.9), that corrects the Faddeev-Popov distribution cut-off at the 
Gribov horizon for over-counting inside the Gribov horizon. 

[We may mention here an alternative approach. The Landau gauge is the singular 
limit a of more regular gauges, and contains a non-local effective drift force KgtefF) eq. 
(4.9). For this reason it may be preferable to calculate with gauge parameter a finite, so 
the drift force, K,, = + a'^Df^d ■ A, remains local, and there is no horizon outside 
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of which the probabihty distribution vanishes exactly. In this case the SD equation (5.5) 
for the effective action F gets replaced by 



The gluon propagator is given by = i?, where R = —j^K[A + V(A)-£^]\a^o, as 

in eq. (5.7). One would hope to solve the SD equations for the full propagators in this 
approach, and not just their infrared asymptotic limit. An advantage of this approach 
is that the solution for finite value of the gauge parameter a could be directly compared 
with the numerical lattice data of [7] and [8] that is taken with stochastic gauge fixing 
and gauge parameter a = 0.1. To control ultraviolet divergences, it will be necessary to 
develop Ward-type identities appropriate to this scheme. They were not needed in the 
present calculation because no ultraviolet divergences appeared in the infrared limit. Such 
identities in BRST-form are available in the 5-dimensional scheme that is based on the 
time-dependent diffusion equation [22], [23], [24], and alternatively one may attempt to 
solve the SD equations of the 5-dimensional scheme non-perturbatively.] 

In the second part of the article, where we calculated the infrared critical exponents, we 
have however ignored the new term K2 in order to compare with other authors, and because 
it is not without interest to calculate the infrared critical exponents non-perturbatively in 
Faddeev-Popov theory with a cut-off at the Gribov horizon. 

It is noteworthy that all our values for the critical exponents in d = 2, 3 and 4 
dimensions agree with exact results for a probability distribution that is cut-off at the 
Gribov horizon namely, the vanishing [17] of the gluon propagator D{k) as k ^ 0, 
and the enhancement [26], [27] and [28] of the ghost propagator [k'^G{k)]~^ 0, (except 
for the first solution in (i = 2 which is marginal, with ac = and qd = —1). The vanishing 
of D{k) at /c = is counter-intuitive, and has no other explanation than the proximity of 
the Gribov horizon in infrared directions. This suppresses the infrared components A{k) of 
the gluon field, and thus of the gluon propagator D{k) = {\A{k)\'^). Since our calculation 
involves a truncation of the SD equations which is an uncontrolled approximation, the 
stability of our results should be tested by estimating corrections. As for the future, an 
immediate challenge is to include the effect of the new term ^2, eq. (4.9), that was not 
evaluated in the present calculation. One must also introduce quarks. 

We wish to compare our values of the infrared asymptotic dimensions with those 
reported by [13], [14], and [15]. But first we must verify whether they also selected the 




(8.1) 
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solution of the SD equations of Faddeev-Popov theory that vanishes outside the Gribov 
horizon. Note that to obtain a particular solution it is suflBcient to require any one of 
its properties, provided that this requirement selects a unique solution. Indeed a unique 
solution was obtained in [13] by requiring that both the gluon and ghost propagators 
D{k) and G{k) be positive. These properties by no means follow from the Faddeev-Popov 
weight (2.1) that oscillates in sign, whereas restriction to the Gribov region does imply the 
positivity of both G{k) and D{k). So in fact the restriction to the Gribov region is also 
implemented in this way in [13]. Likewise the assumptions made in [14] and [15] to obtain 
a solution of the SD equations are equivalent to the horizon condition, eq. (6.4), that we 
imposed in sec. 6. 

It is reassuring that the values given in eq. (7.7) for d = 4 agree qualitatively with the 
values reported in [13], namely ac = [61 — y^l897)]/19 ~ 0.92 and = —2ac ~ —1.84, 
in the sense that the gluon propagator D{k) ~ l/(/c^)^"'""^ vanishes at /c = 0, and the ghost 
propagator G{k) ~ l/(/c^)^+"<3 is enhanced. This may be an indication that these qualita- 
tive features of the solution are not merely an artifact of the approximations made. For the 
two treatments of the SD equations are quite different. Indeed in [13], the gauge condition 
is treated off-shell, by imposing the Slavnov- Taylor identities to determine longitudinal 
parts of vertices, and by using the method of [25] to adjust the gluon propagator. On 
the other hand we have treated the gauge condition on-shell, so only transverse quantities 
occur. There is a similar qualitative agreement for d = 4 with the values reported in [14], 
ac = (77— a/(2281))/38) 0.769479, and ac = — 2aG, where an angular approximation 
was made. The approximations made in [15] appear to be similar to ours, although the 
method of solution is quite different. The value reported there for d = 4, ac = 1 and 
an = —2, agrees with our first solution.^ 

We also wish to compare our results with numerical Monte Carlo studies of propagators 
in Landau gauge. Numerical gauge fixing to the Landau gauge is achieved by minimizing, 
with respect to gauge transformations, the lattice analog of FA{g) — J cfix which 
indubitably produces configurations that lie inside the Gribov horizon. This gauge fixing, 
like stochastic gauge fixing, has a Euclidean weight that is everywhere positive, without 
over-counting. However it is not in the class of Faddeev-Popov gauges for which the 
determinent alternates in sign, so a comparison with analytic calculations by the Faddeev- 
Popov method does not have a completely clear interpretation. 

^ After completion of this article, L. von Smekal has kindly informed me that the first value, 
aa = [93 — y^l201)]/98 in d = 4 dimensions, was also obtained by C. Lerche [29]. 
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The infrared behavior of the lattice propagators is very sensitive to finite-volume 
effects, and control of the volume dependence at fixed /5 = IN/qq is required. In particular 
D{k) does not and should not vanish at A; = at any finite lattice volume, but only 
when extrapolated to infinite volume. We have not attempted here to fit the data of [7] 
and [8] without an estimate of the finite-volume correction and the effect of the finite 
gauge parameter, but this is a promising avenue for future comparison of numerical and 
analytic results. However we do note that it was reported in [7], with stochastic gauge 
fixing at gauge parameter a = 0.1 (with Landau gauge at a = 0), that a fit to the Gribov 
formula, D{k) = -|- M"^]"^, (strong infrared suppression) can explain the gross 

feature of the data. Recent studies in Landau gauge at finite lattice volume indicate a 
suppression of the gluon propagator in the infrared [30], and are not incompatible with 
an enhancement of the ghost propagator [31]. The infrared behavior of the lattice gluon 
propagator D{k) has been studied in SU(2) gauge theory in Landau gauge in d = 3 
Euclidean dimensions [32]. It was found that D{k) has a maximum at 350 MeV 
(normalized to the physical value of the string tension) that is practically /9-independent, 
and that D{k) decreases as k decreases below this value. This decrease is interpreted 
as resulting from the proximity of the Gribov horizon in infrared directions. A similar 
behavior is expected for the 3-dimensionally transverse part of the gluon propagator in 
Coulomb gauge, in 4 Euclidean dimensions. This has been observed, and an extrapolation 
to infinite lattice volume at fixed /3 was in fact found, notably, to be consistent with the 
vanishing of D(/c) at /c = [33] and [34]. We emphasize that this behavior is not seen at 
finite lattice volume but only in the extrapolation to infinite lattice volume, at fixed (3. 
For this reason it is important to extend the lattice calculations in Landau, Coulomb and 
stochastic gauges to larger volumes, and to extrapolate to infinite lattice volume before 
attempting a fit to continuum formulas. 

So what have we learned about propagators and the confinement problem in QCD? 
We may summarize results qualitatively by the statement that in the infrared region in 
non-perturbative Landau gauge there is strong suppression or vanishing of the would-be 
physical gluon propagator, and strong enhancement of unphysical propagators.^*^ This is 

The "unphysical propagator" that is infrared-enhanced may be either the ghost propagator 
in non-perturbative Landau-gauge Faddeev- Popov theory, or the 44-component of the gluon prop- 
agator in Coulomb gauge. The ghost propagator coincides, approximately, with the remnant of 
the longitudinal gluon propagator that survives the Landau-gauge limit in stochastic quantization. 
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true both for the analytic solutions of the Schwinger-Dyson equations obtained in [13], 
[16], [15], [29] and here, and for the numerical lattice data just discussed, with similar 
numerical data in Coulomb gauge. We expect that these qualitative features will stand 
the test of time. They provide a simple intuitive picture of confinement in which the sup- 
pressed massless physical gluon disappears from the physical spectrum while the enhanced 
unphysical components provide a long-range color-confining force. (This long-range force 
should also confine quarks, but that has not been addressed here.) As discussed previously 
[35] both features may be understood as the result of the restriction to the Gribov re- 
gion, that results from the identification of gauge-equivalent configurations. The infrared 
suppression of the transverse gluon propagator results from the proximity of the Gribov 
horizon in infrared directions, while the enhancement of the unphysical components is an 
entropy effect that results from high population in the neighborhood of the Gribov horizon, 
where the inverse Faddeev-Popov operator is enhanced. 
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Appendix A. Evaluation of integrals 

To evaluate the gluon self-energy Id(p-g)i eq- (6.14), we write 

/•oo 

l/[(p- fc)Y+'*° = r-^(l + aG) / cZx a;"°exp[-a;(j9- fc)2], (A.l) 

Jo 

and similarly for l/[{pf'Y'^'^'^ . This gives 

/•oo nOO 

In{aG) = N[{d-l){k^)^^+^-^T\l + aG)]-' dx dy {xy^^J (A.2) 

Jo Jo 



where 



J = {2n)-'^ J d'^p \p^k^ - {p ■ kf] exp[-a;(p - kf - yp^] 
= id-l) e [2 (47r)^/2 {x + yY+'^/Y^exp[-{x + y)-^xyk'^] 



(A.3) 



We introduce the identity 1 = J dj d{x + y — and change variable according to x — jx' 
and y — jy'. This gives, after dropping primes, 

loiao) = N{e)-^^+''^\2 (477)^^/2 ^2^^ ^ aoT^K (A.4) 
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nOO pOO pOO 

K= dx dy d-f 5{x + y ^ 1) {xy)'^° -f^'^°-'^/^ exp[--fxyk^] 
Jo Jo Jo 

^ (.2.-2aa-i+d/2 r{2aG + l-d/2) T^-aG + d/2) 

^ ' T{d - 2aG) 

This gives 

N V{2aG + l-d/2)V^{-aG + d/2) 
~ 2 (47r)d/2 r2(l + aG)T{d - 2aG) 

where we used ckd = —2aG — (4 — (i)/2. 

To evaluate the ghost self-energy, Ig{olg)i eq. (6.15), we use the identities 



1 



[(/c)2]2+ai, r(2 + 



1 /"°° 

-/ dx x'^'^+^eM-xk^), 

^ocd)Jo 



111 f 

X dz exp{-y[z{k - pf + (1 - z)k^]}, 
Jo 

which allow us to cancel the leading power of k explicitly. This gives 



where 



This gives 



(p2)i+«G T{2 + aD) r(l + aG) 

X / dx dy dz x'^''+^ L 
Jo Jo Jo 

L = (27r)-'^ J d'^k \p^k^ - {p ■ kf] (/ -2p-k) 

X e^p{—[x + y(l — z)]k^—yz{k — p)^}. 

^ (p^)^ - 1) x + y-2yz / yz[x + y{l - z)]p'^ \ 
2 {A-K)d/^ {x + y)^+d/^ ^""P V x + y A 



iV(rf-l) (p2)l-^C 



Tdx Hdy dz X--+' y--+' ; 

Jo Jo Jo (X + ?/)2+'i/2 

^ / _ y^[a^ + y(l -^)]f 

V x-\-y 
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We again introduce the identity 1 = J S{x + y — ^) and change variables according to 
X = ^x' and y = 7?/'. This gives, after dropping primes, 

POO poo POO pi 

J= d-f dx dy dz S(x + y-l) a;'*^+^ (1 - 2yz) 

Jo Jo Jo Jo (A. 14) 

X 7"^*° exp ^ - yz[x + y{l - z)]^p'^^ , 

J ^ (p2)aG-ir(i _ f dy f dz (1 - (1 - 2yz) 

Jo Jo (A. 15) 

X [yz{l-yz)r^-\ 

where we again used ao = — 2q!g — (4 — (i)/2. We change variable of integration to 
u = yz{l — yz), with du = y{l — 2yz)dz, and obtain 

J = (p^o-l ^-1 r(i _ f' dy y^-o (1 _ yyaa-l+d/2 

Jo 

^ (p2)aG-i r(l-aG) r(l + 2aG) T{-aG + d/2) 
aa r{aG + 1 + d/2) 

where we have again used ao = —2aG — (4 — d)/2. This gives 

N{d-1) TT 



(A.16) 



Ig{o.g) = 



2 (47r)'^/2 sin(7rQ;G) 

r(2aG + l) V{-aG + d/2) 



(A.17) 



X 



r2(aG + 1) r(-2aG + d/2) T{aG + 1 + d/2) ' 
where we have used T[aG) r(l — ckg) = tt/ sin(7rQ;G)- This integral is positive for ckg < <^/4. 
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